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Preface 


Carbon  dioxide  lasers  have  become  the  topic  of  increased 
research  over  the  last  few  years.  One  area  of  active  inves- 
tivation  is  the  plasma  chemistry  occurring  in  the  laser 
discharge.  These  studies  are  ultimately  related  to  the 
development  of  flowing  and  nonflowing  sealed  C02  lasers. 

Such  laser  systems  offer  versatility  which  is  advantageous 
to  the  Air  Force  and  industry.  It  is  hoped  that  the  present 
study  aids  in  the  understanding  of  plasma  chemistry  pro- 
cesses in  the  CC>2  laser  discharge. 

I would  like  to  thank  the  whole  staff  of  the  Advanced 
Concepts  Group  of  the  Air  Force  Aero  Propulsion  Laboratory. 
Special  thanks  goes  to  my  advisor,  Alan  Garscadden,  for 
his  guidance  and  Michael  Stamm  for  the  chemistry  code  used 
in  the  theoretical  predictions.  Thanks  is  also  due  to 
Richard  Franklin  for  his  helpful  suggestions  on  the  code 
and  draft  of  this  report.  Last  but  not  least,  I would 
like  to  thank  my  wife,  Hattie,  for  her  patience  and  for 
preparing  most  of  the  figures  used  in  this  report. 
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Abstract 


The  neutral  chemistry  of  the  following  mixtures  have 
been  studied  by  using  a quadrupole  mass  spectrometer  and 
a theoretical  model:  C02/Ar,  N2/Ar,  C02/N2/Ar,  C02/He/Ar, 
and  C02/N2/He/Ar.  C02  dissociation  occurs  by  direct  elec- 

tron impact  and  appears  to  depend  on  initial  C02  concentra- 
tion and  mixture.  Very  little  N2  dissociation  occurs  in 
a pure  N2  discharge;  however,  the  dissociation  of  N2 
increases  when  it  is  added  to  a mixture  containing  C02. 

This  dissociation  is  less  than  1%.  Consequently,  the 
dominant  nitrogen  oxide  in  mixtures  containing  CC>2  and 
N2  is  NO,  and  the  concentration  of  NO  never  exceeds  1%. 

The  species  CHO  and  CH20  are  present  at  high  discharge 
currents  if  there  is  a small  amount  of  water  contamination 
in  mixtures  containing  C02« 
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e,  SYNTHESIS  OF  THE  PLASMA  CHEMISTRY 

OCCURRING  IN  HIGH  POWER  C02  LASERS 

| I.  Introduction 

[ 

jj  Background 

I The  most  important  problem  facing  closed-cycle  C02 

[ lasers  is  that  under  sealed-off  conditions,  the  original 

\ C02/N2/He  mixture  changes  considerably  due  to  the  plasma 

[ chemistry  occurring  in  the  laser  discharge.  Reactions 

of  prime  concern  are  the  dissociation  of  C02  into  CO  and 
O,  the  formation  of  nitrogen  oxides,  and  the  formation 
of  significant  amounts  of  negative  ion  species.  The  irre- 
versible loss  of  C02  reduces  the  output  power,  while  the 
nitrogen  oxides  and  negative  ion  impurities  alter  the 
properties  of  the  laser  discharge. 

Two  methods  are  used  to  solve  the  problem  of  C02  dis- 
sociation (Ref  12:2039).  One  method  is  to  employ  heated 
platinum  or  copper  as  a catalyst.  The  other  method  is 
to  use  the  catalytic  action  of  H20  or  H2.  Small  amounts 

of  water  vapor  or  H2  in  the  discharge  reform  the  C02 

through  the  following  reactions: 

CO  + OH  + C02  + H (1.1) 

CO  + H02-*-  C02  + OH  (1.2) 
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The  OH  comes  from  the  dissociation  of  H20,  and  the  H02 
comes  from  the  dissociation  of  H2  ar.^  subsequent  combina- 
tion of  H with  02. 

The  nitrogen  oxides  that  are  formed  quench  the  C02 
vibrational  excitation  and  also  compete  with  the  C02  for 
vibrational  excitation  by  electron  inelastic  collisions. 

These  impurities  also  affect  the  laser  gain  and  discharge 
impedance  (Ref  4).  The  negative  ion  population  can  become 
comparable  to  (and  in  some  cases  even  exceed)  the  electron 
population/ resulting  in  discharge  arcing  (Ref  13:268). 

The  effects  of  the  nitrogen  oxides  and  negative  ion  species 
mentioned  here  limit  long-term  operation  of  closed-cycle 
systems. 

Previous  experimental  studies  of  the  plasma  chemistry 
have  usually  involved  pure  C02  (Refs  1,  7,  and  8)  or  the 
total  laser  mixture  consisting  of  C02/N2/He  (Refs  2,  11, 

12,  18,  22,  and  24).  C02  dissociation  in  a pure  CC>2  dis- 

charge does  not  necessarily  occur  at  the  same  rate  as  C02 
dissociation  in  a mixture.  Although  it  is  a simple  job 
to  track  the  buildup  of  certain  species  in  a C02/N2/He 
discharge,  the  identification  of  the  primary  reactions 
involved  is  a complicated  task.  In  order  to  get  a clear 
picture  of  the  reactions  involved,  it  is  advantageous  to 
start  initially  with  simple  mixtures.  Then,  known  amounts 
of  additional  gases  are  titrated  in.  The  above  procedure 
represents  the  "synthesis"  approach  in  studying  the  chem- 
istry of  the  laser  discharge.  Important  initial  steps 


toward  this  type  of  study  have  been  made  by  Buser  and  Sul- 
livan (Ref  5) , and  Smith  and  Austin  (Ref  20) . 

Theoretical  models  of  the  chemistry  have  also  been 
made.  For  example,  Shields  and  Smith  (Ref  18)  and  Shields 
et  al.  (Ref  19)  have  modeled  the  negative  ion  effects  in 
C02  convection  lasers  and  TEA  C02  lasers,  respectively. 

While  these  models  provide  some  insight  into  the  chemistry 
occurring  in  the  discharge,  they  still  rely  on  important 
assumptions  regarding  reaction  rates.  In  many  of  the  models 
for  example,  the  reverse  reaction  rates  of  most  reactions 
were  set  equal  to  zero. 

Review  of  Previous  Studies 

The  experimental  and  theoretical  studies  of  the  plasma 
chemistry,  some  of  which  were  mentioned  in  the  last  sec- 
tion, will  now  be  discussed  in  detail.  The  results  of 
each  study  will  be  given,  and  if  applicable,  a comparison 
with  similar  studies  will  be  made. 

In  March  1969,  Corvin  and  Corrigan  (Ref  7)  investi- 
gated the  dissociation  of  C02  by  electron  impact  in  the 
positive  column  of  a glow  discharge.  It  was  reported  that 
from  the  measurement  of  the  rate  coefficients,  the  dissoci- 
ation takes  place  via  electron  impact.  The  conclusion 
was  that  the  dominant  reaction  is ; 

C02  + e"  -*•  CO  + 0 + e~  (1.3) 

Four  months  later,  Gasilevich  et  al.  (Ref  11)  reported  that 
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equilibrium  for  reaction  (1.3)  is  reached  in  several  seconds 
for  large  diameter  discharge  tubes  and  fractions  of  seconds 
in  small  ones.  They  went  on  to  say  that  C02  dissociation 
becomes  greater  as  pressure  and  flow  velocity  are  decreased 
and  discharge  current  is  increased. 

The  following  February,  Buser  and  Sullivan  (Ref  5) 
studied  initial  processes  in  C02  glow  discharges.  They 
reported  that  not  only  does  collisional  dissociation  con- 
tribute to  the  C02  dissociation  process  as  in  reaction 
(1.3),  but  dissociative  attachment  as  well.  This  dissoci- 
ative attachment  reaction  is  given  by 

C02  + e"  CO  + 0“  (1.4) 

Buser  and  Sullivan  also  reported,  for  the  first  time,  that 
H2  in  the  mixtures  reduces  the  dissociation  process. 

The  experimental  studies  described  so  far,  except 
the  one  by  Corvin  and  Corrigan,  involved  a flowing  gas 
system.  In  April  1970,  Karube  and  Yamaka  (Ref  12)  investi- 
gated C02  dissociation  in  a non-flowing  sealed  CC>2  laser. 
Dissociation  of  C02  molecules  was  found  to  be  consistent 
with  the  reaction  , 

C02  + CO  + -|  02  (1.5) 

with  an  equilibrium  constant  determined  by  the  discharge 
parameters.  For  example,  there  was  a 60%  dissociation 
for  a discharge  current  of  10  mA,  and  the  dissociation 
increased  slightly  with  increasing  current. 

4 
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Austin  and  Smith  (Ref  2)  identified  the  major  posi- 
tive ions  in  the  total  C02  laser  mix  and  pure  C02.  The 
cathode  region  was  sampled  with  a time-of-f light  mass 
spectrometer,  and  the  major  ions  for  the  total  mix  were 

C02+,  °2+/  co+'  N2+'  an<*  H20+*  The  Pure  C02  discharge 
had  the  following  ions:  C02+,  0 2+,  C0+,  H20+,  0+,  and 
C+.  Although  these  results  are  valid  in  the  negative 
glow  and  cathode  fall  regions,  one  should  note  that 
they  may  not  necessarily  be  true  in  the  positive  column 
of  the  glow  discharge  where  most  of  the  laser  excitation 
occurs.  Austin  and  Smith  published  the  results  of  their 
study  in  1972. 

In  January  1973,  Wiegand  and  Nighan  (Ref  24)  used 
a theoretical  rate  equation  model  and  reported  that  numer- 
ous secondary  minority  species,  such  as  those  composed 

of  carbon,  oxygen,  and  nitrogen  are  produced  in  signifi- 

-4 

cant  quantities  on  a time  scale  of  10  sec.  The  dissoci- 
ation products  CO  and  0 are  produced  on  the  same  time 
scale.  These  neutral  minority  species  can  reach  concentra- 
tions in  excess  of  the  charged  particle  densities,  and 
they  are  found  to  be  directly  related  to  the  dominant 
ionic  species.  For  typical  laser  discharge  parameters, 
the  negative  ion  species  can  become  comparable  to  the 
electron  concentration.  The  dominant  negative  ion  species 
are  C03",  C04”,  N02”,  0~,  and  N03“.  Shields  and  Smith 
(Ref  18)  also  reported  five  years  later  that  in  a typical 
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COj/Nj/He/I^O  mixture,  the  principal  ions  are  C03  , CO^ 
and  H-.  These  negative  ions  are  greater  than  ten  percent 
of  the  electron  concentration,  but  never  exceed  it.  How- 
ever, if  the  gases  are  recycled  or  if  there  is  an  air  leak, 

the  dominant  negative  ions  are  NO  ~ and  NO--.  These  ions 

2 ^ 

can  exceed  the  electron  concentration.  The  results  men- 
tioned here  are  only  applicable  to  C02  convection  lasers, 
in  TEA  C02  lasers,  Shields  et  al.  (1976)  stated  that  in 
Sealed  C02/N2/He/02/C0/H2/H20  mixtures,  the  dominant  ions 
are  C04~  and  H-^  However,  in  systems  with  little  02  or 
HjO  added,  the  dominant  ion  is  C03_  (Ref  19) . The  studies 
cited  in  this  paragraph  used  a computer  model  of  the  plasma 
chemistry.  Prince  and  Garscadden  (1975)  experimentally 
determined  that  N0^~ , N03  , 02~,  0-,  and  C03  are  indeed 
the  principal  negative  ions  formed  (Ref  15) . However,  there 
is  some  discrepancy  concerning  the  dominant  negative  ion. 
Prince  and  Garscadden  found  that  the  dominant  negative 
ion  is  N02  while  the  theoretical  models  predict  CC>3  . 

They  did  find  that  the  C03~  concentration  rises  sharply 
at  low  discharge  currents,  and  at  higher  currents,  N02 
becomes  the  dominant  negative  ion.  i 

Smith  and  Austin  (Ref  20)  continued  the  study  of  C02 
dissociation  in  pulsed  and  continuous  lasers  in  August 
1973.  This  study  measured  C02  dissociation  for  a variety 
of  E/p  values.  Important  steps  were  made  in  studying  C02 
dissociation  for  a variety  of  mixtures  with  C02  as  a 


; 
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component.  The  dissociation  is  greater  for  C02/N2/He  than 
for  C02/N2»  and  the  dissociation  is  greater  for  C02/N2 
than  for  pure  C02. 

Thoenes  and  Kurzius  expanded  the  research  from  C02 
dissociation  to  the  study  of  the  plasma  chemistry  important 
in  the  formation  of  impurities  (Ref  22) . This  study,  done 
in  August  1976,  concentrated  on  closed-cycle  C02  lasers. 

The  two  mixtures  studied  were  He/N2/C02  and  H2/N2/C02. 

The  study  by  Thoenes  and  Kurzius  was  again  theoretical. 

Two  studies  in  November  1976  dealt  mainly  with  posi- 
tive and  negative  ion  reactions  in  carbon  dioxide,  and 
ion-molecule  reactions  in  C02  and  C02/C0  mixtures.  These 
were  done  by  Alger  and  Rees  (Ref  1) , and  Coxon  and  Moruzzi 
(Ref  8) , respectively.  Alger  and  Rees  reported  that  for 
pure  CO 2,  the  dominant  ions  are  02+  and  C02+.  They  also 
measured  the  reaction  rate  for  the  following  charge  exchange 
reaction  involving  these  ions: 

°2+  + CO,  -*■  C02+  + O,  (1.6) 

Coxon  and  Moruzzi  measured  the  rates  for  the  formation 
of  positive  C02  cluster  ions  over  a wide  range  of  E/n  values 
The  C02  cluster  ion  (C02.C02)  is  formed  by  the  following 
reactions : 

C02+  + 2C02  j (C02.C02)+  + C02  (1.7) 

Other  than  the  studies  by  Buser  and  Sullivan  (Ref  5) ? 
and  Smith  and  Austin  (Ref  20) , there  has  been  no  important 


plasma  synthesis  study.  This  type  of  study  will  almost 
surely  be  advantageous  in  understanding  plasma  chemistry 
processes  in  CO2  laser  discharges. 

Problem  Statement  and  Approach 

The  present  problem  is  to  obtain  an  accurate  picture 
of  the  plasma  chemistry  occurring  in  a C02  laser  discharge 
The  "synthesis"  approach,  as  described  previously,  will 
be  used  in  order  to  reduce  the  ambiguity  of  the  results. 
The  neutral  chemistry  of  the  following  mixtures  will  be 
studied  experimentally  by  using  a quadrupole  mass  analyzer 
C02/Ar,  N2/Ar,  CC>2/N2/Ar,  C02/He/Ar,  and  C02/N2/He/Ar. 

The  volume  flow  ratio  of  components  in  the  multi-component 
mixtures  above  will  reflect  the  same  volume  flow  ratio 
found  in  typical  C02  lasers  (i.e.,  as  in  a 1:1:8  mixture 
of  C02/N2/He) . A theoretical  model  of  each  mixture  will 
be  developed  and  compared  to  the  experimental  data. 
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II.  Experimental  Arrangement  and  Procedure 


The  complete  mass  analyzer  system,  not  to  scale,  is 
shown  in  Pig.  1.  The  test  gases  were  flowed  through  a 
discharge  tube  connected  on-line  to  a quadrupole  mass 
spectrometer.  In  the  data  acquisition  system,  there  was 
no  provision  for  automatic  reduction  of  the  data.  The 
data  reduction  had  to  be  done  by  hand  calculation  and  with 
the  aid  of  a computer.  It  is  the  purpose  of  this  section 
to  explain  each  of  the  parts  of  the  total  system  shown 
in  Fig.  1.  The  data  reduction  process  is  outlined  in  Appen- 
dix A. 

Gas  Handling  System 

A,  B,  C,  and  D in  Fig.  1 represent  the  gas  bottles 
where  the  four  gases  of  this  experiment  were  stored.  These 
gases  were  CO 2,  He,  N2,  and  Ar;  and  they  were  rated  at 
a purity  of  99.998%.  The  regulated  flows  were  measured 
by  Hastings  mass  flowmeters.  The  measured  flows  were  accu- 
rate to  ±1%  for  He  and  Ar,  and  ±5%  for  N2  and  C02.  The 
entry  of  the  gases  into  mixing  manifold (E)  was  controlled 
by  two  valves:  an  on-off  valve  and  a needle  valve.  The 
needle  valve  was  used  to  control  the  flow  rate  of  gases 
into  the  mixing  manifold.  The  needle  valves  held  much 
of  the  back  pressure  from  the  gas  regulators,  and  this 
pressure  was  kept  between  14.7  and  18  psi  in  order  to 
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Experimental  Arrangement.  For  explanation,  see  text. 


avoid  damage  to  these  valves. 


The  mixing  manifold  consisted  of  a hollow  cylinder. 
Normally,  a packing  of  glass  beads  is  used  to  aid  in  a 
thorough  mixing  of  the  gas  components.  In  this  experiment, 
glass  beads  were  not  imperative  because  of  the  slow  flow 
vised.  The  flow  of  each  component  gas  in  a mixture  was 
less  than  10  standard  cubic  centimeters  per  minute  (SCCM) , 
and  the  total  flow  was  around  12  SCCM.  At  these  slow  flow 
rates,  the  gases  have  sufficient  time  in  the  glass  tubing 
leading  to  the  discharge  tube  (F  in  Fig.  1)  to  mix  before 
entering  the  discharge  tube  itself  (G  in  Fig.  1) . The 
mixing  was  also  aided  by  the  low  pressure  used  in  the  exper- 
ment  (2  torr) . This  is  not  the  case  for  very  high  flow 
rates.  Another  advantage  of  a slow  flow  is  that  the  total 
flow  rate  of  a mixture  of  gases  is  just  the  sum  of  the 
flow  rates  of  the  individual  components. 

The  Discharge  Tube 

The  discharge  tube  (G)  was  of  pyrex  construction  with 
an  outer  jacket  for  cooling  water.  The  inner  diameter 
was  2.50  cm,  and  the  length  was  approximately  30-40  cm. 

The  hollow  cylindrical  electrodes  (H)  , which  were  separated 
by  a distance  of  30  cm,  were  made  of  Kovar. 

The  discharge  was  controlled  by  a current  regulated 
high  voltage  power  supply  (I) , and  the  cathode  of  the  dis- 
charge tube  was  located  upstream  of  the  gas  flow.  The 
maximum  current  attainable  was  around  80  mA.  In  a Nj/Ar 
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discharge,  arcing  occurred  at  about  70  mA.  The  sampling 
orifice  (J)  was  located  midway  between  the  two  electrodes 
so  that  it  sampled  from  the  positive  column  of  the  glow 
discharge.  The  orifice  used  in  this  experiment  was  2 mils 
in  diameter.  The  residence  time  (see  Appendix  B)  is  the 
length  of  time  that  a portion  of  the  gas  stays  in  the  dis- 
charge before  it  is  sampled  by  the  spectrometer.  It  is 
given  by 


where  A is  the  cross-sectional  area  of  the  discharge  tube, 

1 is  the  length  from  the  cathode  to  the  sampling  orifice, 
p is  the  pressure  of  the  gas,  pa  is  the  standard  pressure 
(760  torr) , and  H is  the  total  flow  rate  given  in  SCCM. 

The  residence  time  for  this  experiment  was  chosen  to  be 
1 sec.  This  choice  permitted  a flow  that  was  measurable 
by  the  flowmeters.  Once  a residence  time  is  chosen,  a 
restriction  is  made  on  the  total  flow  rate  H for  a par- 
ticular gas  pressure  p.  Solving  for  H in  equation  (2.1) 

and  using  the  parameters  that  are  applicable  to  this  experi 
. 2 

ment  (i.e.,  Jl=15cm,  A=ira  where  a is  the  inner  radis  of 
the  discharge  tube,  and  T=lsec) , the  flow  rate  is 

H = (5 . 81) p (2.2) 

The  pressure  in  the  discharge  was  measured  by  a dif- 


i 


ferential  capacitance  manometer  (K)  manufactured  by  MKS 


r 

Instruments,  Inc.  All  pressure  measurements  were  made 
upstream  of  the  gas  flow  behind  the  cathode  of  the  dis- 
[ charge  tube.  The  reference  pressure  was  obtained  via  a 

i 

[ cold  trap  (L)  at  liquid  nitrogen  temperature.  The  procedure 

for  obtaining  the  reference  pressure  is  as  follows.  First, 
the  manometer  head  and  cold  trap  are  pumped  down  by  the 
same  forepump  that  evacuates  the  discharge  tube.  To  do 
this,  valves  VI  and  V2  are  opened.  Next,  valve  VI  is  closed 
and  the  cold  trap  is  filled  with  liquid  nitrogen.  Now 
the  manometer  is  at  the  pressure  of  the  cold  trap,  and 
the  meter  is  zeroed.  To  measure  the  pressure  of  the  dis- 

I 

charge  tube  relative  to  the  cold  trap  pressure,  valve  V2 
is  closed  and  VI  opened.  It  is  important  to  keep  the  cold 
trap  filled  with  liquid  nitrogen.  If  the  trap  is  allowed 
to  increase  in  temperature,  then  the  reference  pressure 
changes  which  results  in  a change  in  measurement  of  the 
source  pressure. 

The  Quadrupole  Mass  Spectrometer 

The  spectrometer  used  in  this  experiment  was  a quadru- 
pole mass  filter  manufactured  by  Extranuclear  Laboratories, 
Inc.  The  model  is  an  EMBA  II.  The  molecular  species  to 
be  analyzed  are  first  ionized  before  detection.  The  com- 
plete mass  spectrometer  system  is  designated  M in  Fig  1. 
There  are  two  differential  pumping  stages.  The  first  is 
pumped  by  a turbomolecular  pump,  and  the  second  by  a dif- 
fusion pump.  The  first  stage  is  usually  pumped  down  to 
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a pressure  of  2.4  x 10  torr  while  the  second  stage  to 

— 8 

a pressure  of  8.4  x 10  torr.  The  molecular  beam  from 
the  orifice  is  chopped  by  an  electric  tuning  fork  chopper 
before  entering  the  rear  orifice  (N) . This  chopped  beam 
is  then  ionized  by  the  ionizer  (0) , and  the  resultant  ion 
beam  is  focused  by  a series  of  lenses  (P)  before  entering 
the  quadrupole  region  (Q) . 

The  lens  system  consists  of  a lens  to  control  the 
ion  energy,  an  extractor,  and  three  focusing  lenses.  The 
first  lens  is  given  a positive  potential,  and  only  positive 
ions  of  a certain  minimal  kinetic  energy  can  surmount  the 
potential  barrier.  The  extractor  is  given  a negative  poten- 
tial to  accelerate  the  ions  toward  the  focusing  lenses. 

The  three  focusing  lenses  are  given  potentials  (positive 
or  negative)  to  focus  the  ion  beam  into  the  quadrupole 
region  (Q) . 


Theory  of  Quadrupole  Mass  Filtering.  The  most  com- 
plete theory  of  quadrupole  mass  filtering  is  given  by  1 
et  al.  (Ref  14).  However,  the  following  summary  is  adapted 

from  Ref  16.  Four  symmetric  hyperbolic  electrodes,’ 

2 2 

described  by  the  equation  X -Y  = ±C  where  C is  a constant, 
make  up  the  quadrupole  mass  filter.  Two  opposite  elec- 
trodes are  given  the  potential  ({> (t)  = U+Vcos (wt) , where 
0 is  a dc  voltage  and  V is  the  peak  amplitude  voltage  of 
a radio-frequency  source  of  angular  frequency  w.  The  other 


14 


sign  of  the  first  two.  Fig.  2 shows  this  configuration. 
Under  these  conditions,  a quadrupole  field  is  set  up  with 
the  field  on  axis  equal  to  zero. 


If  an  icr  is  injected  with  a velocity  parallel  to 
the  z axis,  it  undergoes  a transverse  motion.  The  equa- 
tions of  motion  in  the  xy  plane  are  given  by  Mathieu  equa- 
tions, and  the  solution  is  either  exponential  or  oscil- 
latory. For  a given  U and  V,  ions  of  a particular  q/m 
will  have  stable  trajectories  through  the  quadrupole 
assembly.  This  stable  trajectory  consists  of  oscillatory 
motions  about  the  z axis  coupled  with  a velocity  parallel 
to  the  z axis  such  that  the  particle  does  not  strike  the 
electrodes.  Since  the  pressure  in  the  back  chamber  of 

_7 

the  spectrometer  is  less  than  10  torr,  the  mean  free 
path  of  the  ions  is  great  enough  so  that  collisions  among 
ions  are  negligible  when  undergoing  the  oscillatory  motions. 
All  ions  without  the  right  q/m  value  will  ultimately  strike 
the  electrodes  and  be  removed  from  the  ion  beam.  The  mass 
resolution  is  given  by 


_m_  _ 0.126 

Am  “ 0. 16784-U/V 


» (2.3) 


Note  that  the  resolution  of  the  mass  filter  can  be  varied 
by  varying  the  ratio  U/V.  If  U/V  = 0.16784,  then  a 
theoretical  infinite  mass  resolution  is  achieved.  For 
infinite  mass  resolution. 
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Figure  2.  Quadrupole  Mass  Filter 


U - 1.212  mf2rQ2  (2.4) 

V = 7.219  mf2rQ2  (2.5) 

where  m is  the  mass  in  amu,  f is  the  frequency  of  the  rf 

source  in  MHz,  and  r is  the  radius  of  the  circle  inscribed 

o 

tangent  to  the  four  hyperbolic  electrodes  (see  Fig.  2)  . 

In  order  for  mass  separation  to  occur,  the  ion  of 
the  incorrect  q/m  must  stay  in  the  quadrupole  long  enough  so 
that  it  can  be  filtered  out.  Hence,  this  requirement  puts 
a restriction  on  the  ion  energy  and  length  L of  the  elec- 
trodes. The  approximate  relationship  is  given  by 


p i «2_  ^ / Ann  ^ 

Eion  < 25  f L ( m )m 


(2.6) 
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where  E.  _ is  the  kinetic  energy  of  the  ion  in  eV,  and 
L is  the  length  of  the  quadrupole  assembly  given  in  cm. 
Lighter  ions  will  be  transmitted  more  than  heavier  ions 
because  of  their  greater  velocity.  Although  lighter  and 
heavier  ions  may  have  the  same  kinetic  energy,  their  respec- 
tive velocities  through  the  quadrupole  are  different.  There- 
fore, there  is  a transmission  discrimination  against  heavier 
ions.  To  compensate  for  this,  the  ion  energy  control  may 
be  increased  with  increasing  mass  units.  The  spectrometer 
used  in  this  experiment  did  not  have  the  provision  of 
increasing  ion  energy  with  increasing  mass  units. 

There  is  also  a maximum  aperture  radius  of  ion  injec- 
tion into  the  quadrupole  region.  If  the  ion  is  injected 
too  far  off  from  the  z axis,  it  will  strike  the  elctrodes 
even  if  it  has  the  correct  q/m.  For  ions  injected  parallel 
to  the  z axis,  the  maximum  aperture  radius  of  ion  injection 
to  insure  that  none  will  strike  the  electrodes  is  given 
by 

a * 4 r (2.7) 

3 o m 


From  equations  (2. 3) , (2.6)  and  (2.7);  it  is  evident 
that  the  resolution  and  transmission  of  the  ions  through 
the  filter  can  be  controlled  by  varying  the  applied  volt- 
ages D and  V.  The  achievement  of  good  resolution  and 
transmission  depends  on  the  stability  of  these  two  voltages. 
The  quadrupole  mass  filter  used  in  this  experiment  has 
provisions  for  adjusting  the  resolution  and  transmission 
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by  varying  the  ratio  of  the  two  voltages.  Also,  it  can 
be  seen  from  equations  (2.4)  and  (2.5)  that  mass  selection 
can  be  done  by  either  varying  U and  V,  or  f.  In  the  Extra- 
nuclear  quadrupole  system,  U and  V are  varied  while  f is 
held  constant.  Therefore,  the  two  voltages,  U and  V,  play 
an  important  role  in  ion  transmission  and  mass  discrimina- 
tion. 


Particle  Detector.  All  Extranuclear  particle  detectors 
are  designed  so  that  the  first  dynode  is  not  on  axis  with 
the  mass  filter.  The  particle  detector  is  denoted  by  R 
in  Fig.  1.  This  arrangement  prevents  photons  and  neutrals 
from  striking  the  first  dynode  producing  random  noise  in 
the  electronics.  The  paraxial  detector,  in  which  the  multi- 
plier is  offset  but  parallel  to  the  axis  of  the  quadrupole, 
is  used  in  the  EMBA  II  mass  spectrometer. 

After  exiting  the  mass  filter,  the  ions  enter  a space 
in  which  the  surrounding  walls  are  at  ground  potential 
(i.e.,  the  spectrometer  case  is  grounded).  The  partition 
that  separates  the  mass  filter  from  the  particle  detector 
contains  a circular  aperture  which  approximately  matches 
the  viewing  area  of  the  first  dynode  of  the  particle  multi- 
plier. The  first  dynode  is  usually  kept  at  a high  negative 
voltage  so  that  the  ions  are  drawn  to  it.  The  multiplier 
voltage  was  kept  at  -4  kVin  this  experiment. 

Quadrupole  Controls . The  quadrupole  control  unit 
(S  in  Fig.  1)  serves  several  functions.  Circuitry  is 
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provided  to  control  the  voltage  output  of  U and  V to  a 
fixed  value.  U and  V may  also  be  swept  over  certain  ranges 
with  the  aid  of  a sweep  generator;  and  the  sweep  is  over 
a specified  range  of  m/q.  The  output  of  the  sweep  was 
connected  to  the  x input  of  the  x-y  plotter.  There  are 
also  provisions  for  indicating  the  value  of  m/q  that  the 
quadrupole  is  tuned  to.  The  resolution  mode  can  be  varied 
by  letting  U=yV-6  in  equation  (2.3)  where  y and  6 are  con- 
stants. The  constants  y and  6 can  be  controlled,  and  these 
are  the  resolution  and  Am  controls,  respectively . The  resolu- 
tion was  set  at  a ratio  of  8.15,  and  the  Am  was  set  at 
+1.8  volts.  These  settings  eliminated  peak  splitting  in 
the  output  sweep.  Another  control  to  improve  peak  shaping 
is  called  the  pole  bias;  and  this  control  was  set  at  -0.4 
volts  which  improved  the  symmetry  of  the  peak.  The  resolu- 
tion, Am,  and  pole  bias  controls  depended  somewhat  on  the 
size  of  the  signal  and  sensitivity  setting. 

Detection  Electronics 

The  detection  electronics  consisted  of  a particle 
multiplier,  lock-in  amplifier,  and  an  x-y  plotter.  These 
parts  are  schematically  represented  by  R,  T,  and  U in  Fig. 

1.  The  particle  multiplier  has  been  discussed  in  some 
detail  in  a previous  section.  As  stated  earlier,  a tuning 
fork  chopper  is  located  in  front  of  the  second  aperture 
to  the  back  chamber  of  the  spectrometer.  This  chopper 
chopped  the  molecular  beam,  thus  providing  a modulated 


signal  to  the  particle  detector.  The  signal  from  the  par- 
ticle detector  was  fed  to  the  lock-in  amplifier  with  the 
frequency  of  the  tuning  fork  as  a reference.  The  signal 
output  from  the  lock-in  amplifier  was  then  led  to  the  y 
input  of  the  x-y  plotter. 

The  settings  on  the  amplifier  and  x-y  plotter  were 
kept  the  same  on  all  spectra  with  the  exception  of  an 
increase  in  amplifier  sensitivity  on  some.  The  time  con- 
stant of  the  amplifier  was  300  ms,  and  the  sensitivity  was 
5 mV  on  the  following  spectra:  C02/Ar,  N2/Ar,  and  C02/N2/Ar. 
The  spectra  of  those  mixtures  which  contained  helium  were 
taken  at  a sensitivity  of  2 mV.  The  2 mV  sensitivity  was 
used  on  the  following  mixtures:  C02/He/Ar  and  C02/N2/He/Ar . 
This  higher  sensitivity  produced  a recognizable  helium 
signal.  The  period  of  the  sweep  on  all  recorded  spectra 
was  300  sec. 

Special  Calibrations 

Aside  from  the  instrumentation  calibration  given  in 
the  EMBA  II  manual,  there  are  two  important  non- instrument 
calibrations  that  must  be  made.  These  are  the  contact 
potential  offset  of  the  ionizer  energy  meter  and  the  spectro 
meter  function.  Both  of  these  calibration  factors  play 
an  important  role  in  the  data  reduction  process  (see  Appen- 
dix A) . These  factors  will  be  discussed  in  some  detail 
below. 

Contact  Potential.  The  energy  reading  from  the  meter 
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is  not  the  energy  of  the  ionizing  electrons.  The  meter 
reading  is  just  the  potential  difference  between  the  fila- 
ment and  grid,  but  the  energy  of  the  electrons  is  actually 
distributed  about  this  value.  Aside  from  this  fact,  there 
is  an  electron  energy  offset  between  the  ionizer  and  the 
meter  due  to  the  connections  between  the  two.  It  is  a 
well  known  fact  that  there  is  a contact  potential  differ- 
ence between  two  dissimilar  metals  connected  together. 

This  potential  difference  is  expected  to  be  on  the  order 
of  a few  volts. 

To  determine  this  contact  potential,  appearance  poten- 
tial plots  are  run  for  different  gases.  The  difference 
between  the  reading  from  the  plot  and  the  actual  appearance 
potential  is  the  amount  of  offset.  In  this  experiment, 
appearance  potential  plots  were  run  for  C02,  N2,  He,  and 
Ar.  The  offset  was  approximately  3.5  volts.  This  voltage 
offset  also  takes  into  account  the  effect  of  the  space 
charge  of  the  ionizer  and  the  energy  loss  in  overcoming 
the  work  function  of  the  cathode. 

Spectrometer  Function.  The  mass  spectrometer  used 
in  this  experiment  is  a number  density  measuring  device. 

As  stated  earlier,  the  test  gas  is  first  ionized  before 
detection.  Fig.  3 shows  the  construction  of  the  ionizer. 

The  sample  molecular  beam  is  led  into  the  cylindrical  grid 
where  it  is  crossed  by  an  electron  beam.  It  is  constructive 
at  this  point  to  derive  an  expression  for  the  total  ion 
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Figure  3.  Electron  Impact  Ionizer 


current  density  that  is  formed  by  the  ionizer.  Out  of 
this  expression,  the  spectrometer  function  will  be  defined. 
The  simplified  ionizer  shown  in  Fig.  4 will  be  considered 
in  the  following  analysis. 

The  number  of  ionizing  collisions  per  unit  volume 
per  unit  time  between  electrons  and  molecules  is  given 


ft  = n n < Vo  (V) > 
e s s 


(2.8) 


where 


« the  electron  number  density 

= the  molecular  number  density  of  molecular 
species  s 

= the  magnitude  of  the  relative  velocity 
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Figure  4.  Simplified  Ionxzer 


o_(V)  * the  total  ionization  cross-section  of  molecular 
species  s 

The  <>  denotes  a calculated  average  in  velocity  space  of 

the  product  Vo  (V) . 

s 

It  is  assumed  that  the  electrons  possess  most  of  the 
kinetic  energy  such  that  V ~ Ve  where  is  the  speed  of 
the  electrons.  It  is  further  assumed  that  the  speed  of 
the  electrons  is  constant.  Therefore,  equation  (2.8) 
reduces  to 


n n a V 

esse 


(2.9) 
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The  magnitude  of  the  interacting  electron  current  density 
in  the  ionizer  is  assumed  to  be  given  by 


*eneve 


(2.10) 


where  q is  the  magnitude  of  the  electronic  charge.  That 
is,  the  assumption  has  been  made  that  the  ionization  pro- 
cess is  not  having  to  compete  with  other  processes  such 
as  inelastic  scattering,  dissociation,  etc.  Solving  for 
nfi  in  equation  (2.10)  and  substituting  into  equation  (2.9), 
the  following  result  is  obtained  for  the  number  of  ions 
per  unit  volume  per  unit  time  formed  by  the  ionizer: 


n = 


ns°s 


(2.11) 


The  number  of  ions  formed  per  second  in  the  ionizer  is 
obtained  by  multiplying  equation  (2.11)  by  the  volume  of 
the  interaction  region: 

» - 7 D2X,i  - | D2L  ns°s  12-121 

where  D is  the  diameter  of  the  entrance  aperture  and  L is 
the  length  of  the  ionizer  (see  Fig.  4) . If  these  ions 
are  now  swept  out  as  fast  as  they  are  produced,  then  the 
ion  current  density  becomes 

«2-»» 

Equation  (2.13)  states  that  the  extracted  ion  current  density 
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is  proportional  to  the  product  of  the  molecular  number 
density  and  total  ionization  cross-section  of  molecular 
species  s.  Equation  (2.13)  also  assumes  that  the  molecules 
are  singly  ionized. 

When  the  quadrupole  is  tuned  to  a particular  mass 
number,  the  signal  is  proportional  to  equation  (2.13). 

(PH)S  - (Cs jeL> o#n#  <2.14) 

where  (PH) s stands  for  the  signal  measured  in  terms  of 

the  deflection  of  the  chart  recorder,  and  C is  a constant 

s 

of  proportionality.  C is  in  general  different  for  each 
molecular  species.  Equation  (2.14)  defines  what  is  known 
as  the  spectrometer  factor.  The  spectrometer  factor  is 
given  by 


Fs  = CsjeL  (2.15) 

Using  the  definition  of  the  spectrometer  factor,  equation 
(2.14)  becomes 


(PH)  = Fan  (2.16) 

s s s s 

The  factors  in  equation  (2.15)  describe  certain  charac- 
teristics of  the  spectrometer.  These  factors  will  be  dis- 
cussed in  some  detail  below. 

The  electron  current  density  (j  ) and  the  length  of 

© 

the  interaction  region  (L)  are  both  characteristics  of 
the  ionizer.  The  electron  current  density  of  the  ionizer 
depends  on  the  rate  of  thermionic  emission  from  the  filament, 
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and  it  also  depends  on  the  potential  difference  between 
the  filament  and  the  grid.  The  interaction  length  is  just 
the  length  of  the  ionizer. 

Cg  includes  the  characteristics  of  the  particle  multi- 
plier and  detection  electronics.  The  response  of  the 
particle  multiplier  is  a function  of  the  number  density 
and  impact  energy  of  the  atoms  or  molecules  of  a particular 
species.  In  the  effusive  flow  regime,  the  velocity  distri- 
bution of  the  molecules  is  Maxwellian.  Although  a flowing 
gas  mixture  is  used  in  this  experiment,  it  is  assumed  that 
the  velocity  distribution  of  each  species  is  still  Max- 
wellian about  a coordinate  frame  moving  with  the  same 
velocity  as  that  of  the  gas.  The  response  of  the  particle 
multiplier  will  be  the  same  for  all  species  with  the  same 
mass  number  if  equal  number  densities  of  these  species 
strike  the  particle  multiplier.  In  the  regime  of  a free- 
jet  expansion,  the  molecules  reach  a terminal  velocity 
at  a point  of  about  30D  where  D is  the  diameter  of  the 
orifice  (Ref  22:12).  This  terminal  velocity  is  given  by 

V = ( (2y/y-1)  (kT  /m))**  ' (2.17) 

s 

where  m is  the  molecular  mass,  Tg  is  the  temperature  of 
the  source,  k is  Boltzmann's  constant,  and  y is  the  ratio 
of  the  specific  heat  at  constant  pressure  to  that  at  constant 
volume.  If  the  flow  contains  a mixture  of  gases,  then 
V is  the  mass  average  velocity  and  m is  the  average  mass. 

The  ratio  of  the  specific  heats,  y,  is  just  that  of  the 
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mixture.  Therefore,  molecular  species  in  a mixture  have 
the  same  terminal  velocity.  Again  the  signal  is  the  same 
for  all  molecular  species  with  the  same  mass  number  (if 
their  number  densities  are  equal) . It  should  be  noted 
that  while  two  molecular  species  may  have  the  same  mass 
number,  their  internal  degrees  of  freedom  may  be  different. 
The  number  of  secondary  electrons  emitted  by  the  first 
dynode  of  the  particle  multiplier  depends  on  the  total 
energy  which  is  the  sum  of  the  translational  and  internal 
energies.  Therefore,  it  is  assumed  that  the  internal  energy 
is  negligible  compared  to  the  translational  energy  in  order 
to  simplify  the  data  reduction  process.  This  is  a valid 
assumption  since  ions  are  accelerated  by  the  high  negative 
voltage  of  the  first  dynode. 

There  are  advantages  and  disadvantages  to  both  of 
the  flow  regimes  described  above.  While  a correct  picture 
is  given  of  the  chemical  makeup  of  the  source  in  the  effu- 
sive flow  regime,  the  sensitivity  of  the  detection  elec- 
tronics must  be  high  in  order  to  obtain  recognizable  sig- 
nals. This  high  sensitivity  is  required  due  to  the  weak 
intensity  of  the  molecular  beam.  The  free-jet  expansion 
is  advantageous  in  producing  a high  intensity  beam,  but 
it  has  the  disadvantage  of  altering  the  original  chemistry 
of  the  source.  The  reason  is  that  a free-jet  expansion 
is  collision  dominated,  and  chemical  reactions  are  bound 
to  occur  within  the  beam  before  detection.  The  most  desir- 
able flow  regime  is  neither  effusive  nor  free-jet  expansion. 
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any  predictions  concerning  transition  flow  are  limited 
to  the  analytical  models  of  effusive  flow  and  free-jet 


expansion.  There  is  little  knowledge  to  date  concerning 

the  regime  of  transition  flow. 

The  number  density,  n , in  the  ionizing  region  of 

the  ionizer  obviously  depends  on  the  partial  pressure  of 

molecular  species  s in  the  source.  If  the  molecular  beam 

consists  of  a mixture  of  gases,  then  diffusive  separation 

can  occur  (Ref  17) . The  factor  C also  takes  into  account 

s 

this  diffusive  separation  of  components  in  the  beam.  More- 
over, it  will  be  assumed  that  the  fractional  amount  of 
the  total  number  density  of  each  component  in  the  beam 
is  the  same  as  that  in  the  source. 

When  F is  solved  for  in  equation  (2.16) , the  follow- 
s 

ing  result  is  obtained: 

<PH) 

F_  = (2.17) 

*sns 

A plot  of  F vs.  atomic  mass  unit  is  known  as  the  spectro- 
meter function.  To  generate  points  for  this  plot,  differ- 
ent gas  mixtures  are  flowed  individually  through  the  dis- 
charge tube  with  the  number  density  of  each  molecular 
species  known.  From  the  experimentally  determined  signal, 
and  the  known  ionization  cross-section  corresponding  to 
the  energy  of  the  ionizer,  an  F corresponding  to  the  atomic 
mass  unit  of  each  molecular  species  can  be  found.  The 
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, accuracy  of  the  plot  is  increased  if  several  individual 

gases  are  used. 

Fig.  5 shows  a plot  of  F vs.  atomic  mass  unit.  The 

following  gases  were  used  to  obtain  this  plot  by  using 

the  method  described  in  the  previous  paragraphs  Hfi,  N2, 

Ar , and  C02.  The  electron  energy  of  the  ionizer  was  set 

at  70  eV.  In  equation  (2.11),  (PH)  is  in  units  of  cm 

(i.e.,  the  deflection  in  cm  of  the  chart  recorder)  and 

2 

of  is  in  atomic  units  (ira  , where  a is  the  Bohr  radius)  . 
s o o 

A close  examination  of  equation  (2.17)  indicates  that 
F is  highly  dependent  on  the  partial  pressure  of  molecular 
species  s in  the  source.  The  reason  for  this  dependency 
is  that  the  number  density  in  the  molecular  beam  increases 
*4  with  an  increase  in  partial  pressure  and  vice  versa.  Under 

a discharge  current  of  80  mA,  the  pressure  increases  by 
as  much  as  0.4  torr  in  C02/Ar.  F is  assumed  not  to  change 
much  for  this  small  increase  in  pressure  due  to  the  dis- 
charge. 

From  the  arguments  given  in  the  Spectrometer  Charac- 
teristics section  concerning  the  constant  of  proportional- 
ity C , it  is  not  unreasonable  to  assume  that  F is  the 

s 

same  for  all  atomic  or  molecular  species  with  the  same 
or  nearly  the  same  mass  number.  This  assumption  allows 
an  F to  be  read  for  any  additional  molecular  species  pro- 
duced in  the  discharge  from  a predetermined  plot.  Because 
of  the  large  dispersion  shown  in  Fig.  5,  F should  be  cal- 
culated separately  for  each  species  in  the  cold  gas  mixture 
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ATOMIC  MASS  UNIT'S 

Figure  5.  Spectrometer  Function 


and  used  in  the  reduction  process.  The  curve  in  Fig.  5 
should  only  be  used  to  obtain  an  F for  the  additional 
molecular  species  produced  in  the  discharge. 


III.  Experimental  Results 


The  experimental  results  of  the  mixtures  considered  in 
this  experiment  will  be  giver)  in  this  section.  In  all 
of  the  mixtures,  the  amount  of  Ar  was  10%  in  order  to  moni- 
tor the  density  changes  during  the  discharge.  The  voltage 
and  pressure  effects  of  the  discharge  are  first  discussed. 
Afterwards,  the  neutral  chemistry  of  the  mixtures  will 
be  presented.  Where  possible,  a comparison  will  be  made 
with  other  similar  experiments  in  the  literature. 

Voltage  Effects 

Fig.  6 is  a plot  of  discharge  voltage  vs.  current 
for  all  of  the  mixtures.  The  sustaining  voltage  of  the 
discharge  depends  on  many  factors.  As  evident  from  Fig. 

6,  it  depends  on  the  particular  gas  mixture.  Other  fac- 
tors include  ionization  cross  sections,  the  reduced  field 
(E/n) , and  the  ionization  efficiency.  Arcing  occurred  in 
the  N2/Ar  discharge  above  a discharge  current  of  approxi- 
mately 75  mA.  One  reason  for  the  high  voltage  of  the  N2/Ar 
discharge  is  that  the  ionization  efficiency  of  N2  is  low 
compared  to  the  other  gases  used  in  this  experiment. 

Pressure  Effects 

The  pressure  in  the  discharge  tube  was  2 torr  for 
all  the  mixtures  before  the  discharge  was  turned  on.  After 
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Figure  6.  Discharge  Voltage  Vs.  Current  of 
Experimental  Mixtures 


the  discharge  was  turned  on,  the  pressure  increased  in 
COj/Ar , C02/N2/Ar,  COj/He/Ar,  and  C£>2/N2/He/Ar;  but  decreased 
by  one-one  hundreth  of  a torr  in  the  case  for  N2/Ar.  This 
pressure  decrease  has  also  been  observed  by  Karube  and 
Yamaka  (Ref  12:2036).  The  pressure  increase  is  due  to 
the  formation  of  new  species  in  the  discharge  and/or 
gas  heating. 

Fig.  7 shows  the  effect  of  the  discharge  on  the  source 
pressure.  The  curves  for  C02/Ar,  COj/Nj/Ar,  C02/He/Ar, 
and  C02/N2/He/Ar  exhibit  the  same  shape  as  the  experimental 
curves  indicating  the  rise  in  dissociation  products  as 
a function  of  discharge  current.  The  buildup  of  certain 
species  as  a function  of  discharge  current  will  be  shown 


Figure  7.  Pressure  Vs.  Discharge  Current  of 
Experimental  Mixtures 

later  in  this  section.  Therefore,  it  is  concluded  that 
the  pressure  effects  were  predominantly  due  to  the  formation 
of  secondary  species. 

Further  support  is  given  to  this  conclusion  by  noting 
the  curve  for  Nj/Ar  in  Fig.  7.  Experimentally  determined 
in  this  thesis,  there  is  little  dissociation  in  a N2/Ar 
discharge.  If  the  conclusion  that  the  rise  in  source  pres- 
sure was  due  primarily  to  the  formation  of  secondary  species 
is  valid,  then  the  curve  for  N2/Ar  is  entirely  consistent 
with  this  conclusion.  Moreover,  the  pressure  rise  in  N2/Ar 
should  be  greater  than  the  other  mixtures  if  gas  heating 
was  the  cause  (i.e.,  the  N2/Ar  discharge  is  very  hot  and 
the  pressure  should  rise  dramatically) . 
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90%  CO2/10%  Ar 

Fig.  8 summarizes  the  chemistry  of  a CO^/Ar  discharge. 
Initially  with  a C02  concentration  of  90%,  the  amount  of 
CO 2 dropped  rapidly  with  discharge  current  until  a current 
of  about  50  mA  was  reached.  The  number  of  C02  molecules 
dissociated  was  greater  in  a C02/Ar  discharge  than  in  any 
other  mixture  considered  in  this  experiment.  The  fractional 
dissociation  was  up  to  44%.  Buser  and  Sullivan  have  shown 
that  dissociation  up  to  80%  is  common  (Ref  5:479). 

This  large  dissociation  can  be  explained  by  the  fol- 
lowing. Smith  and  Austin  reported  in  1974  that  dissociation 
is  preceded  by  the  production  of  an  excited  state  of  C02 
by  electron  impact  (Ref  20:321): 

C02  + e~  -*■  C02(3it)  + e"  (3.1) 

Subsequent  dissociation  occurs  though  the  following  reaction 

C02(3tt)  + e"  -*■  CO  + 0 + e“  (3.2) 

Shields  et  al.  (Ref  19:1589)  gave  the  overall  rate  con- 

-9  3 -1 

stant  for  this  dissociation  mechanism  as  1 x 10  cm  s 
(300°K) . This  rate  is  considerably  larger  than  the  rate 
for  dissociative  attachment: 

C02  + e"  + CO  + 0"  (3.3) 

-13  3 -1 

which  has  a rate  constant  at  300K°  of  5 x 10  cm  s 
(Ref  19:1589).  In  a predominantly  CO,  discharge,  there  are 


Figure  8.  Number  Density  Vs.  Discharge  Current  i 

of  the  Major  Species  in  a 90%  CO2/10%  Ar  : 

Discharge  (Experimental)  1 
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more  interactions  between  electrons  and  CO2  molecules  than 

mixtures  containing  CO2  in  which  the  other  components  are 

comparable  to  the  C02  percentage.  Thus,  there  are  more  C02 

3 

molecules  in  the  ir  state  and  subsequent  dissociations 
through  reaction  (3.2) . 

At  a current  of  80  mA,  the  CO  percentage  was  about 
55%.  This  is  consistent  with  the  44%  dissociation  of  the 
C02  as  noted  earlier.  The  rise  in  CO  percentage  was  con- 
siderably higher  in  C02/Ar  than  in  the  other  mixtures. 

The  apparent  reason  is  the  same  as  that  given  above. 

The  dispersion  of  the  experimental  CO  data  points 
is  inferred  from  the  dispersion  of  the  spectrometer  factor 
given  in  Fig.  5.  Since  the  C02  was  an  initial  component 
in  all  experimental  mixtures,  a spectrometer  factor  was 
obtained  for  it  in  the  cold  gas  and  used  in  the  reduction 
of  the  C02  percentage.  Hence,  dispersion  in  the  C02  per- 
centage cannot  be  determined.  Dispersion  in  the  02  per- 
centage cannot  be  determined  because  the  spectrometer  fac- 
tor for  02  was  not  taken  from  actual  experimental  data, 
but  from  the  curve  in  Fig.  5.  Therefore,  the  dispersion 
in  percentage  can  only  be  given  in  the  case  of  CO. 

The  only  explanation  offered  for  the  decrease  in  02 
at  high  currents  is  that  at  a current  of  50  mA,  the  dis- 
sociation of  O2  occurs.  With  the  dissociation  of  02»  the 
percentage  of  atomic  oxygen  should  increase.  Atomic  oxygen 
should  be  readily  detectable  (Ref  5:474);  however,  none 
was  detected  in  the  C02/Ar  discharge.  Furthermore,  0 was 
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not  detected  in  either  of  the  mixtures  studied  in  this 


experiment.  Two  mechanisms  may  contribute  to  the  loss 
of  atomic  oxygen.  First,  the  recombination  of  atomic  oxygen 
on  the  walls  of  the  discharge  tube  to  form  02  is  occurring. 
The  rate  constant  for  such  a process  depends  on  the  diffu- 
sion of  0 to  the  walls.  Secondly,  the  adsorption  of  0 
to  the  walls  and  sampling  orifice  of  the  mass  spectrometer 
is  also  a possible  mechanism  by  which  0 is  lost.  However, 
not  all  of  the  atomic  oxygen  is  lost  by  the  two  processes 
above.  The  remaining  0 in  the  discharge  that  did  not 
recombine  or  get  adsorbed  may  be  below  the  sensitivity 
of  the  spectrometer. 

Ozone  may  also  be  present,  but  was  not  detected. 
Increasing  the  sensitivity  of  the  spectrometer  to  observe 
0 and  03  increases  the  background  noise  of  the  detection 
electronics  causing  the  0 and  0^  signals  to  be  indistinguish- 
able. To  prove  that  0 and  are  present  in  the  discharge, 
spectroscopic  techniques  must  be  employed. 

90%  N2/10%  Ar 

There  was  no  detectable  dissociation  of  N2  in  a N2/Ar 
discharge.  This  conclusion  is  based  on  the  fact  that  there 
was  no  rise  in  the  N signal  and  no  change  in  the  N2  signal. 
This  conclusion  is  also  supported  by  the  pressure  effects  of 
the  N2  discharge  as  shown  in  Fig.  7.  An  increase  in  pres- 
sure would  indicate  the  formation  of  additional  species, 
but  no  such  increase  is  evident  from  the  N2/Ar  curve  in 
Fig.  7. 
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Since  there  was  no  measurable  dissociation,  it  is 
interesting  to  know  where  most  of  the  power  (possessed 

mainly  by  the  electrons)  from  the  discharge  was  transferred. 

— 16  2 

Below  an  E/n  of  10  x 10  Vcm  , most  of  the  power  is  trans- 

— 16  2 

f erred  to  vibrational  excitation.  Above  10  x 10  Vcm  , 
most  is  transferred  to  exciting  the  electronic  metastable 
state  of  N2  and  total  electronic  excitation.  The  plot  for 
fractional  power  transferred  in  N2  vs.  E/n  is  given  in  Fig.  9. 

The  electric  field  in  the  positive  column  of  a N2/Ar 
glow  discharge  can  be  calculated  approximately  in  the  fol- 
lowing way.  By  knowing  the  material  composition  of  the 
electrodes,  the  cathode  fall  of  potential  in  volts  can  be 
determined.  The  composition  of  the  electrodes  was  Kovar, 
and  for  a N2  discharge  the  cathode  fall  is  on  the  order  of 
200  V (Ref  23:229).  From  Fig  6,  the  maximum  and  minimum 
tube  voltages  were  1300  and  1250  V respectively.  The  differ- 
ence between  the  tube  voltage  and  the  cathode  fall  gives 
the  voltage  V across  the  positive  column.  The  electric  field 
of  the  positive  column  is  then  V/l  where  1 is  the  length  of 
the  column.  1 will  be  taken  to  be  the  length  between  the 
anode  and  cathode.  Hence,  i 


1100V  _ „ -1 

= 36.7  Vcm 


max  30  cm 


E_,_  = = 35  Vcm-1 


min  30  cm 


(3.4) 

(3.5) 


As  a first  approximation,  the  gas  density  will  be 
taken  to  be  the  density  without  a discharge.  At  a pressure 
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Figure  9.  Fractional  Power  Transfer  in  a N»  Discharge  (Ref  3) 


f 


1 

of  2 torr  and  a temperature  of  293°K,  the  density  of  n=2.505 
X 1019  (■—)  = 6.59  X 1016  cm"3  (see  Appendix  B)  . Using 
this  value  for  n and  equations  (3.4)  and  (3.5),  the  maxi- 
mum and  minimum  E/n  values  are  given  respectively  by 

E /n  = 5.57  x 10"16  Vcm2  (3.6) 

max 

E . /n  = 5.31  x 10"16  Vcm2  (3.7) 

mm 

At  these  two  values  of  E/n,  over  90%  of  the  power 
goes  into  vibrational  excitation  of  which  some  smaller  frac- 
tion goes  eventually  into  translational  energy.  For  an 
— 16  2 

E/n  of  18  x 10  Vcm  , ionization  accounts  for  only  5 % 
of  the  power  input.  The  rest  of  the  power  is  transferred 
to  the  excitation  of  the  metastable  state  and  total  elec- 
tronic excitation. 

40%  CO2/50%  N2/10%  Ar 

The  chemistry  of  a C02/N2/Ar  discharge  is  illustrated 
in  Fig.  10.  The  initial  C02  concentration  of  40%,  dropped 
to  30%  at  a discharge  current  of  10  mA.  The  maximum  dis- 
sociation occurred  at  a current  of  80  mA  when  the  concentra- 
tion reached  approximately  15%.  The  saturation  point  also 
appeared  to  occur  around  80  mA.  The  fractional  dissociation 
was  62.5%  compared  to  44%  in  C02/Ar  discharge.  Although 
a greater  number  of  molecules  were  dissociated  in  a CO,,  /Ar 
discharge,  the  fraction  of  C02  dissociated  was  greater  in 
the  C02/N2/Ar  mixture.  In  addition  to  C02  dissociation 
by  the  mechanism  illustrated  in  (3.1)  and  (3.2),  there  is 
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possible  dissociation  by  the  collision  of  excited  N2 
molecules.  This  explains  the  greater  fraction  of  C02 
molecules  dissociated  in  the  mixture.  However,  C02  disso- 
ciation by  excited  N2  molecules  is  a hypothesis  which 
requires  rigorous  testing. 

In  a mixture  of  33%  C02  and  67%  N2,  Smith  and  Austin 
(Ref  20:315)  reported  a C02  dissociation  of  20%  above  a 
C02  partial  pressure  of  1.2  torr  and  at  a discharge  current 
of  60  mA.  The  partial  pressure  of  C02  in  the  C02/N2/Ar 
mixture  used  in  this  experiment  was 


pe  = (4.648/11.62)  (2)  = 0.8  torr  (3.8) 

s 

where  Hg  is  the  flow  rate  of  C02,  H is  the  total  flow  rate 
of  the  mixture,  and  p is  the  total  pressure  of  the  mixture 
(see  Appendix  B) . At  60  mA,  the  fractional  dissociation 
was  60%.  This  is  higher  than  Smith  and  Austin's  value  of 
25%.  No  exact  comparison  can  be  made  since  Smith  and  Austin 
gave  no  mention  of  gas  residence  time. 

The  concentration  of  CO  molecules  rose  as  much  as  8.5% 
at  a discharge  current  of  10  mA  and  33%  at  80  mA.  As  in 
the  case  for  C02/Ar,  there  was  no  indication  of  CO  dissocia- 
tion at  high  discharge  currents.  This  observation  was  con- 
firmed when  no  rise  in  the  C signal  occurred.  The  CO  per- 
centage in  the  range  from  10  - 40  mA  is  slightly  higher 
than  that  shown  in  Fig.  10.  The  reason  is  that  CO  and  H2 
have  approximately  the  same  mass  number;  and  in  determining 
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the  CO  signal,  a constant  signal  was  subtracted  from 
the  total  signal  at  mass  number  28.  But  there  was  1%  N2 
dissociation,  and  there  is  no  apparent  reason  why  N2  should 
dissociate  in  a mixture  containing  C02  but  not  in  a N2/Ar 
discharge.  This  same  behavior  was  noticed  in  the  C02/N2/ 
He/Ar  discharge  except  that  the  N concentration  began  to 
increase  slightly  past  40  mA. 

With  this  small  amount  of  N2  dissociation,  various 

nitrogen  oxides  should  be  detectable.  N20  could  not  be 

detected  because  it  has  the  same  or  nearly  the  same  mass 

number  as  C02,  and  the  resolution  of  the  quadrupole  is  not 

sufficient  in  resolving  the  N20  and  C02  signals.  Thus, 

N20  may  have  been  formed.  N02  was  not  formed  in  detectable 

quantities.  This  conclusion  was  based  on  the  fact  that 

a signal  appeared  at  mass  number  46  with  no  discharge.  The 

14 

signal  was  attributed  to  C 02-  If  N02  was  formed,  then 

14 

its  signal  would  increase  the  C 02  signal.  However,  this 
was  not  the  case.  There  was  a trace  amount  of  NO,  and  great 
care  was  taken  in  identifying  it.  In  a C02/Ar  discharge, 
a signal  at  mass  number  30  appeared  only  at  high  discharge 
currents.  Since  there  was  a little  water  contamination 
(<1%) , this  species  was  tentatively  identified  as  CH20. 

It  is  safe  to  assume  that  CH20  appeared  at  high  discharge 
currents  in  the  C02/N2/Ar  mixture  so  that  its  signal  was 
added  to  the  NO  signal.  At  low  currents,  the  signal  at 
mass  number  30  appeared;  and  it  was  attributed  to  NO.  The 
concentration  of  NO  stayed  at  a constant  0.29%  of  the  total 
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ion  current  from  the  spectrometer  in  the  range  of  10  - 


80  mA.  NO  was  not  detected  in  a N2/Ar  discharge  since  there 
was  no  detectable  N2  dissociation. 

The  02  concentration  rose  from  5%  at  10  mA  to  16.5%  at 
60  mA.  Past  60  mA,  the  concentration  appeared  to  decrease; 
however,  more  data  points  are  required  to  reach  a definite 
conclusion.  If  the  02  concentration  indeed  decreased,  then 
the  decrease  was  at  a much  slower  rate  than  that  shown  in 
Pig.  8.  As  can  be  seen  in  Fig.  8,  the  02  concentration 
decreased  rapidly  as  the  discharge  current  was  increased 
past  50  mA  in  a C02/Ar  discharge.  There  are  two  possible 

v 

mechanisms  for  the  slower  decrease  in  a C02/N2/Ar  discharge. 
One  is  that  the  reaction  rate  of  02  formation  on  the  walls 
of  the  discharge  tube  increased  (i.e.,  the  diffusion  of 
0 to  the  walls  increased) , and  the  other  is  that  less  dis- 
sociation occurred  through  electron  impact. 

10%  CO2/80%  He/10%  Ar 

Fig.  11  shows  the  discharge  chemistry  of  a C02/He/Ar 
mixture.  The  purpose  of  the  He  in  C02  lasers  is  to  relax 
the  lower  laser  level  of  the  C02  molecule  to  the  ground 
state.  Thus,  the  He  is  not  expected  to  influence  the 

I chemistry  of  the  discharge  much  (except  in  the  case  of  three 

body  reactions  involving  He) . The  saturation  point  of  C02 
I dissociation  was  approximately  2.5%  concentration  at  a 

f 

| discharge  current  of  90  mA.  This  point  was  determined  by 

[ 1 extrapolating  the  curve  in  Fig.  11  to  high  values  of  the 

t 
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current.  The  fractional  dissociation  was  74%,  and  it  was 


!' 


higher  than  the  dissociation  in  CC^/Ar  and  C02/N2Mr.  At 

40  mA,  the  dissociation  was  65%  which  is  high  compared 

to  the  Gasilevich  et  al.  value  of  60%  dissociation  (Ref 

11:88).  They  had  similar  experimental  conditions  except 

that  the  ratio  of  the  C02  concentration  to  the  He  was 

1:9,  and  the  diameter  of  the  tube  was  22  mm  compared  to 

a diameter  of  25  mm  in  this  experiment.  They  also  reported 

that  the  dissociation  depended  on  the  diameter  of  the  tube 

so  that  this  may  account  for  the  difference  between  their 

values  of  the  dissociation  and  the  value  obtained  in  this 

experiment.  If  C02  dissociation  occurs  through  mechanisms 

(3.1)  and  (3.2),  it  is  difficult  to  explain  why  a greater 

fraction  of  C02  molecules  are  dissociated.  Since  He  relaxes 

the  vibrational  excitation  of  C0o,  there  should  be  more 

*•  I 

C02  molecules  in  the  ground  state.  However,  since  He  has 
a high  ionization  potential  (~  24V)  the  electrons  are 

expected  to  be  more  energetic  in  a predominantly  He  discharge, 
thereby  causing  more  dissociations. 

The  measured  percentage  of  CO  in  the  discharge  was 
higher  than  the  original  C02  concentration  between  35  and 
80  mA  as  shown  in  Fig.  11.  At  10  mA,  the  concentration 
of  CO  reached  6.5%  of  the  original  number  density;  however, 
at  a discharge  current  of  80  mA,  the  concentration  was 
14%.  This  can  be  explained  by  noting  that  the  spectrometer 
function  shows  a large  dispersion  at  mass  28,  and  the  same 
spectrometer  factor  for  CO  was  used  in  the  data  reduction 


miiirr  - 


47 


of  all  the  mixtures  considered  in  this  experiment.  Although 


the  percentages  of  CO  in  Fig.  11  may  be  a little  high, 
the  basic  shape  of  the  curve  is  still  correct.  There  was 
no  problem  with  the  02  concentration  being  higher  than 
the  original  C02  concentration.  The  02  concentration  also 
showed  no  decrease  with  increasing  discharge  current. 

8%  C02/  Hi  N2/  21%  He/  20%  Ar 

The  fractional  dissociation  of  C02  in  a C02/N2/He/Ar 
discharge  was  up  to  76%  as  shown  in  Fig.  12.  The  dissocia- 
tion was  higher  than  any  previous  mixture  considered.  An 
important  trend  was  noted:  as  the  original  C02  concentra- 
tion in  a mixture  was  decreased,  greater  dissociation  of 
CO 2 occurred.  The  partial  pressure  of  C02  was  0.16  torr, 
and  the  dissociation  at  60  mA  was  72.5%.  Smith  and  Austin 
reported  a dissociation  of  68%  for  the  same  conditions 
of  C02  partial  pressure  and  discharge  current  (Ref  20:315). 
No  satisfactory  comparison  between  the  dissociation  reported 
here  and  that  reported  by  Smith  and  Austin  can  be  made 
since  they  did  not  state  their  gas  residence  times,  as 
mentioned  earlier.  , 

One  underlying  mechanism  that  may  explain  the  depen- 
dence of  C02  dissociation  on  C02  concentration  and  mixture 
is  CO  and  O recombination  at  the  walls  of  the  discharge 
tube.  This  recombination  at  the  walls  has  already  been 
introduced  as  one  of  the  reasons  why  atomic  oxygen  was 
undetectable.  Some  of  the  parameters  that  radial  diffusion 


depends  on  are  the  gas  mixture,  the  partial  pressure  of 
the  molecular  species  in  question,  the  total  gas  pressure; 
and  the  temperature.  By  knowing  the  rate  of  CO  and  0 diffu- 
sion to  the  walls  and  the  probability  of  CO  and  0 encounters 
at  the  walls,  a reaction  rate  constant  can  be  computed. 
Recombination  at  the  walls  also  depends  on  tube  geometry 
and  gas  flow.  This  is  a possible  explanation  as  to  why 
Gasilevich  et  al.  found  that  CO2  dissociation  depended 
on  the  tube  geometry  and  gas  flow  (Ref  11) . 

In  Fig.  12,  the  CO  concentration  exceeds  the  initial 
CO  concentration  at  a discharge  current  of  15  mA.  As  noted 
earlier,  this  discrepancy  can  be  attributed  in  part  to 
an  inaccurate  spectrometer  factor  for  CO  in  this  particular 
mixture.  There  was  no  indication  of  saturation  of  the 
CO  percentage  at  high  currents.  However,  small  amounts 
of  carbon  were  detected.  The  percentage  of  C in  the  spectro- 
meter ion  current  fluctuated  from  0.26%  at  10  mA  to  0.087% 
at  80  mA.  The  C may  have  come  from  CO  dissociation  at 
high  currents;  however,  this  is  not  apparent  from  Fig. 

12.  The  dissociation  of  O2  at  high  currents  was  also  not 
detected  as  evident  from  Fig.  12.  The  02  concentration 
saturated  at  8%  of  the  original  number  density  which  was 
the  initial  C02  concentration.  The  high  saturated  concen- 
trations of  CO  and  O2  apparently  are  due  in  part  to  the 
inaccuracy  of  their  respective  spectrometer  factor. 

As  in  the  C02/N2/Ar  mix,  there  was  a small  amount 
of  Nj  dissociation.  The  N concentration  fluctuated  over 
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the  whole  range  of  discharge  currents  similarly  to  the 
C concentration  mentioned  in  the  previous  paragraph.  For 
this  reason,  there  must  have  been  a fluctuation  in  the 
electronics  or  total  gas  flow.  With  this  small  amount 
of  N2  dissociation,  there  was  again  a strong  possibility 
of  nitrogen  oxides  in  the  discharge.  There  was  no  trace 
of  N02,  and  the  presence  of  N20  could  not  be  detected 
because  of  its  similarity  in  mass  number  with  that  of  C02. 

As  stated  in  the  case  of  a C02/N2/Ar  discharge,  the  spectro- 
meter cannot  discriminate  between  N20  and  C02  because  of 
its  limited  resolution.  The  NO  stayed  at  a constant  0.5% 
in  the  range  of  10-80  mA.  This  percentage  was  greater 
than  the  0.2%  reported  in  the  case  of  a C02/N2/Ar  discharge. 
Thus,  the  formation  of  NO  was  greater  in  a C02/N2/He/Ar 
discharge  than  in  a C02/N2/Ar  discharge. 

A signal  at  mass  number  29  was  observed,  and  it  was 
partly  attributed  to  the  formation  of  CHO.  This  CHO  is 
extremely  plausible  since  traces  of  free  carbon  and  atomic 
hydrogen  were  detected.  The  presence  of  atomic  hydrogen 
will  be  explained  in  the  next  section.  The  signal  at  mass 

number  29  was  also  observed  in  the  case  of  a pure  N2  dis- 

14  15 

charge  and  was  tentatively  identified  as  N N . There- 

14  15 

fore,  the  signals  of  CHO  and  N N may  very  well  be  com- 
bined in  the  C02/N2/He/Ar  mixture. 

The  Influence  of  Water  Contamination 

In  all  of  the  mixtures  considered,  there  was  less 
than  one  percent  water  contamination.  The  following 
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additional  species  were  detected  in  all  mixtures  at  low 
discharge  currents:  H20,  OH,  and  H2«  At  high  currents, 
the  following  were  detected:  H20,  OH,  and  H.  The  reason 
that  atomic  hydrogen  appeared  at  high  currents  was  that 
at  these  higher  currents,  dissociation  of  H2  occurred. 

If  water  was  not  present,  then  CHO  and  CH20  in  addition 
to  the  species  mentioned  above  would  not  have  appeared. 
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IV.  Theoretical  Prediction  of  Plasma 
Chemistry  Processes 


In  the  following  theoretical  treatment,  it  will  be 
assumed  that  the  gas  pressure  is  low  enough  so  that  certain 
ternary  reactions  are  negligible.  This  assumption  is  justi- 
fied because  the  gas  pressure  used  in  this  experiment  was 
2 torr,  and  ternary  reactions  at  this  pressure  do  not  con- 
tribute much  to  the  chemistry.  Only  those  with  sizeable 
rates  will  be  considered.  However,  the  bulk  of  the  reac- 
tions considered  in  this  analytical  model  will  be  binary. 


Binary  Reactions 

Consider  the  binary  reaction  m involving  species  Ai 
and  Aj  in  which  species  A^  and  A^  are  formed: 

kf 

Ai  + Aj  f Ax  (4.1) 

kr 

kf  and  kr  are  the  forward  and  reverse  reaction  rates 
respectively.  The  populating  or  depopulating  rate  of  A^ 
is  given  by 

/dlAklX 

m = MVtA.]  - kr[VtAl]  <4’2 

where  the  brackets  denote  number  density.  /dtAj]  \ , 

\ dt  / m 

- and  ~ are  also  equal  to  the  expression 

\ dt  / m'  \ dt  / m 

on  the  right  hand  side  of  equation  (4.2).  Equation  (4.2) 
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can  be  generalized  to  include  the  case  where  more  than 
two  species  are  formed  from  the  reaction  of  with  Aj . 

It  will  be  assumed  that  the  formation  and  depletion 
of  a particular  atomic  or  molecular  species  are  due  entirely 
to  chemical  reactions.  Under  this  assumption,  the  total 
rate  equation  for  species  A^  is 


where  the  sum  is  over  all  reactions  that  form  A.  . The 
result  of  this  formulation  is  a set  of  coupled  nonlinear 
differential  equations  involving  the  number  densities  of 
the  molecular  species  in  the  laser  discharge.  A closed 
analytical  solution  of  these  equations  is  not  feasible. 
Therefore,  they  are  solved  numerically. 

The  treatment  presented  here  on  binary  reactions  can 
easily  be  extended  to  include  ternary  reactions.  The  only 
difference  is  that  since  three  species  are  reacting,  the 
first  term  on  the  right  hand  side  of  equation  (4.2)  includes 
the  product  of  the  three  number  densities.  . . . 

Chemistry  Code 

A list  of  the  important  reactions  and  their  respective 
rate  constants  are  given  in  Appendix  C.  These  were  not 
the  only  ones  considered,  but  they  are  the  ones  that  influ- 
enced the  chemistry  the  most.  A few  of  the  reactions  have 
a temperature  and  average  electron  energy  dependence.  The 
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temperature  was  taken  to  be  300°K,  and  the  average  electron 
energy  was  taken  to  be  1 eV.  Ion  reactions  were  not  con- 
sidered because  it  was  assumed  that  the  neutral  chemistry 
was  perturbed  little  by  them.  Also,  the  reverse  rate  con- 
stants of  the  reactions  in  Appendix  C were  set  equal  to 
zero  because  they  were  not  available.  The  resultant  differ- 
ential equations  were  solved  by  the  modified  Euler  method. 

The  reactions,  their  rates,  and  the  differential  equation 
solver  comprise  the  code  which  was  developed  by  M.  Stamm 
(Ref  21) . 

The  code  required  a value  for  the  electron  number 
density,  and  this  value  was  kept  constant  throughout  the 
temporal  evolution  of  the  species.  To  find  this  value, 
an  approximate  E/n  for  the  mixture  was  calculated  using 
the  same  method  as  that  outlined  in  the  case  for  pure  N2. 
Next,  for  the  value  of  E/n  computed,  the  drift  velocity 
of  the  electrons  was  determined.  Finally,  the  number  density 
of  electrons  was  found  from  the  following  relation: 


ne  qvdA 


(4.4) 


where  I is  the  measured  discharge  current,  A is  the  cross- 
sectional  area  of  the  discharge  tube,  q is  the  magnitude 
of  the  electronic  charge,  and  v^  is  the  drift  velocity 
of  the  electrons.  The  drift  velocity  was  taken  from  the 
curves  in  Appendix  D. 

In  each  mixture,  four  values  of  the  reduced  field 
were  computed  corresponding  to  four  values  of  the  discharge 
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current  (10,  20,  40,  and  80  mA)  . Consequently,  this  gave 
four  values  of  the  electron  number  density.  The  code  was 
run  for  each  value  of  ng,  and  the  species  concentrations 
after  a time  of  1 sec  were  compared  to  the  experimental 
values.  The  1 sec  was  consistent  with  the  experimental 
value  of  the  gas  residence  time. 

It  was  found  that  the  published  value  of  the  rate 
constant  (e.g..  Ref  19)  for  C02  dissociation  by  electron 
impact  did  not  result  in  CC>2  dissociation  in  agreement 
with  experiment.  When  the  published  value  was  used  in 
the  code,  the  C02  concentration  went  essentially  to  zero 
in  1 sec.  The  dissociation  was  far  greater  than  that  experi- 
mentally .observed  . A rate  constant  providing  agreement 
with  the  experiment  was  obtained  in  the  following  way. 

In  the  CC>2  derivative,  the  term  involving  C02  dissociation 
by  electron  impact  is  much  greater  than  any  other  term 
in  the  same  rate  equation.  Hence,  the  CC>2  rate  equation 
can  be  given  approximately  by 


d [CO- ] 

sr kf[co2,ne  . <«-5> 

Integrating  (4.5)  and  applying  the  initial  condition  that 
at  t=0,  [C02)  = tC02)o,  the  following  result  is  obtained: 

{C02J  = IC02]q  e'kfnet  (4.6) 

At  t=l  sec, 

ICO,]  = [CO,l  e"kfne  (4.7) 

2 2 0 
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Solving  for  the  forward  rate  constant  k^,  the  following 
is  obtained: 


[C02]q  and  [CC^]  can  be  obtained  from  experiment  so  that 
an  effective  rate  constant  for  a given  nfi  suitable  to  the 
experiment  can  be  determined.  When  the  rate,  calculated 
from  (4.8),  is  used  in  the  code,  then  the  theoretical  curve 
for  the  CC>2  concentration  vs.  discharge  current  approxi- 
mates the  experimental  curve. 

Theoretical  and  Experimental  Comparison 

As  evident  from  Figs.  8 and  10-16,  there  is  good  agree- 
ment between  the  experimental  and  theoretical  C02  concentra- 
tion. As  mentioned  earlier,  this  agreement  is  due  to  the 
fact  that  the  rate  constant  was  computed  from  the  experi- 
mentally observed  dissociation  assuming  a single  exponen- 
tial decay  and  a constant  electron  number  density.  The 
small  difference  between  the  two  curves  is  due  to  the  fact 
that  the  other  terms  in  the  CC^  derivative  were  neglected. 
The  corresponding  rise  in  CO  shows  fair  agreement  between 
theory  and  experiment  in  all  of  the  mixtures.  It  should 
be  remembered  that  a possible  source  of  error  in  the  experi- 
mental measurement  of  CO  is  the  inaccuracy  of  the  spectro- 
meter factor  used  in  normalizing  the  experimental  CO  per- 
centage. The  calculated  rate  constant  indicates  that  C02 
dissociates  by  the  following  mechanism: 
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Figure  14.  Number  Density  Vs.  Discharge  Current 
of  the  Major  Species  in  a 40%  CO2/50%  N_/10%  Ar 
Discharge  (Theoretical)  " 
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Figure  15.  Number  Density  Vs.  Discharge  Current 
of  the  Major  Species  in  a 10%  CO2/80%  He/10%  Ar 
Discharge  (Theoretical) 
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(4.9) 


C02  + e -►CO  + O + e 


As  indicated  in  Appendix  C,  the  rate  constant  calculated 
from  equation  (4.8)  depends  on  the  mixture  and  is  lower 
than  the  published  value. 

The  model  predicts  that  the  concentration  of  0 exceeds 
(>2  in  all  mixtures.  In  C02/Ar  and  C02/He/Ar,  the  02  is 


lost  by 


02  + e -*•  0 + 0 + e 


(4.10) 


-9  3-1 

with  a rate  constant  of  1 x 10  cm  s (Ref  19) . This 
rate  is  larger  than  the  recombination  rate  of  atomic  oxygen 


given  by 


o + o + m-*-o2  + m 


(4.11) 


—33  6 — i 

where  the  rate  constant  is  7.2  x 10  cm  s (Ref  6), 
and  M stands  for  all  species  in  the  discharge  and  the  wall. 
In  addition  to  reaction  (4.10),  atomic  oxygen  is  formed 


N + 02  -*•  NO  + 0 


N + NO  -*•  N2  + 0 


(4.12) 


(4.13) 


in  the  mixtures  containing  nitrogen.  The  rate  constants 
for  (4.12)  and  (4.13)  are  given  by  8.5  x 10~37  cm3  s 1 
(Ref  6)  and  2.7  x 10-11  cm3  s"1  (Ref  6)  respectively.  The 
reason  that  reactions  (4.12)  and  (4.13)  are  important  is 
due  to  the  large  amount  of  nitrogen  dissociation.  In  a 
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current  range  of  10  - 18  mAf  the  N percentage  rises  from 

11  to  25%  in  C02/N2/Ar  and  5 to  11%  in  CC^/^/He/Ar.  The 

rate  constant  for  N2  dissociation  by  electron  impact  had 

-9  3 -1 

to  be  assumed  (1  x 10  cm  s ) , and  it  may  have  been 
too  high  because  the  high  N concentrations  were  not  observed 
experimentally . 

The  dominant  nitrogen  oxide  in  both  C02/N2/Ar  and 
C02/N2/He/Ar  is  NO  as  observed  experimentally  ( < 1%) . It 
is  formed  by  reaction  (4.12)  and  by 

N + O + M -*  NO  + M (4.14) 

— 32  6 — 1 

where  the  rate  constant  is  1 x 10  cm  s (Ref  6) . In 
spite  of  the  low  rate  constant  of  reaction  (4.14),  the 
importance  lies  in  the  large  concentration  of  N and  O. 

The  additional  reaction  given  by 

0 + N02  -*■  02  + NO  (4.15) 

-12  3 -1 

which  has  a rate  constant  of  9.2  x 10  cm  s (Ref  6) 
is  important  in  C02/N2/Ar.  Since  reaction  (4.15)  has  a 
fast  rate,  N02  does  not  build  up  in  quantities  comparable 
to  NO.  This  is  confirmed  by  the  model.  , 
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V.  Conclusions  and  Recommendations 


COj  dissociation  appears  to  depend  on  initial  C02 
concentration  and  mixture.  Dissociation  is  smaller  with 
greater  amounts  of  initial  C02  and  vice  versa.  One  method 
to  use  in  finding  out  whether  dissociation  depends  on  mix- 
ture is  to  vary  the  initial  concentration  of  C02  within 
one  mixture.  The  ultimate  reason  for  this  behavior  may 
lie  in  the  kinetics  of  the  discharge. 

The  calculated  rate  constant  for  C02  dissociation 
lies  between  those  of  dissociative  attachment  and  dissocia- 
tion by  electron  impact  (Ref  19) . This  calculated  rate 

is  larger  than  that  of  dissociative  attachment  by  a factor 
3 

of  10  . Therefore,  it  is  concluded  that  CC>2  dissociation 
occurs  by  electron  impact.  The  difference  between  the 
published  rate  and  the  calculated  one  can  be  attributed 
to  experimental  error  and  chemistry  changes  in  the  flow  from 
orifice  to  detector.  These  chemistry  changes  undoubtedly 
change  the  concentrations  of  all  reactant  species.  The 
problem  of  chemistry  changes  in  a molecular  beam  can  be 
avoided  by  direct  spectroscopic  measurements  of  the  source 
in  determining  the  concentration  of  a particular  species. 

The  concentration  of  atomic  oxygen  exceeds  that  of 
molecular  oxygen  in  all  mixtures.  The  reason  is  due  to 
the  fact  that  the  dissociation  rate  of  C>2  by  electron  impact 
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is  larger  than  the  recombination  rate  of  0 to  form  02« 

This  conclusion  comes  from  the  model  which  did  not  include 
wall  recombination.  Experimentally,  molecular  oxygen  was 
detected  in  higher  concentrations  than  that  predicted  by  j 

the  model  with  no  traces  of  atomic  oxygen  observed.  The  j 

02  concentration  reaches  a maximum  between  discharge  currents 

of  40  and  60  mA,  and  then  it  starts  to  decrease  with  increas-  j 

1 

ing  current.  This  behavior  of  the  02  concentration,  pre-  j 

dieted  by  the  model,  occurs  in  all  mixtures  except  C02/N2/Ar. 
Experimentally,  the  behavior  of  the  02  concentration  men- 
tioned above  occurred  only  in  the  C02/Ar  discharge. 

Little  N2  dissociation  occurs  in  a pure  N2  discharge. 

Most  of  the  input  power  goes  into  vibrational  excitation 
of  which  a smaller  portion  goes  eventually  into  trans- 
lational energy.  However,  the  dissociation  of  N2  increases 
when  it  is  added  to  mixtures  containing  C07.  This  dis- 
sociation is  less  than  1%,  and  the  rate  constant  is  some- 

-9  3 -1 

what  less  than  the  assumed  value  (1  x 10  cm  s ) . A 
possible  explanation  is  that  in  a mixture  containing  CC>2 
and  N2,  N2  dissociation  occurs  by  the  collision  of  an  excited 
C02  molecule  with  a N2  molecule.  1 

The  dominant  nitrogen  oxide  in  a C02/N2/Ar  or  C02/N2/ 

He/Ar  discharge  is  NO.  In  both  cases,  this  was  observed 
experimentally  and  is  also  predicted  by  the  chemistry  code. 

The  NO  concentration  never  exceeds  1%  in  the  discharge 
current  range  of  10  to  80  mA.  The  concentrations  of  N02 
and  N20  are  essentially  zero  in  this  current  range. 
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If  there  is  less  than  1%  water  contamination,  H20, 

OH,  and  H2  appears  in  all  mixtures  at  low  currents.  Atomic 
hydrogen  is  present  at  higher  currents  due  to  the  dissocia- 
tion of  H2-  In  the  mixtures  containing  C02,  the  additional 
species  CHO  and  CH20  are  present.  The  species  CHO  and 
CH20  can  be  identified  positively  if  the  gas  is  first  dried 
before  entering  the  discharge  region. 

The  kinetic  modeling  can  be  improved  in  two  important 
ways.  First,  a calculation  of  the  electron  energy  distribu- 
tion using  a solution  to  the  Boltzmann  energy  equation 
for  the  E/n  and  mixture  of  gases  to  which  it  pertains  should 
be  made.  Then,  an  integration  of  this  distribution  with 
the  appropriate  inelastic  electron  scattering  cross-sections 
could  give  the  proper  electron  kinetic  rates.  Secondly, 
the  model  should  be  extended  to  include  non-neutral  chemistry 
so  that  ionization  and  recombination  are  properly  treated. 
Such  a model  will  provide  more  accurate  information  on 
the  plasma  chemistry  occurring  in  C02  laser  discharges. 
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.ii. 


Appendix  A 
Data  Reduction 


When  determining  the  true  percentage  of  each  gas  com- 
ponent in  a mixture,  three  factors  must  be  considered. 

These  factors  are  the  spectrometer  characteristics,  the 
density  change  due  to  the  rise  in  gas  temperature  under 
discharge  conditions,  and  the  ionization  cross  section. 

The  spectrometer  function  incorporates  the  spectrometer 
characteristics,  and  it  is  discussed  in  the  experimental 
arrangement  and  procedure  section.  The  remaining  two  factors 
will  be  discussed  below. 

Argon  Correction  Factor 

The  number  density  of  the  source  changes  under  dis- 
charge conditions  due  to  gas  heating.  Since  Ar  does  not 
influence  the  chemistry  of  the  source,  it  can  be  added 

in  small  percentages  to  monitor  this  change.  Ar  may  con-  j 

tribute  to  third  body  reactions,  but  the  pressure  of  the 
source  is  low  enough  such  that  the  probability  of  third 

body  reactions  is  very  small.  The  Ar  correction  factor  ! 

is  defined  by: 

KAr  E (PH,Ar/!,H'Ar  <A'1'  | 

where  (PH)Ar  is  the  "peak  height"  in  cm  deflected  by  the  1 

1 

chart  recorder  of  the  Ar  signal  in  the  cold  gas  (no 

J 
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is  that  same  "peak  height"  in  the 


discharge) , and  (PH) 'Ar 
discharge.  Experiments  have  shown  that  in  general,  K^:^. 
The  reason  being  that  the  number  density  decreases  under 
discharge  conditions.  However,  there  have  been  occasions 
when  KAr  was  less  than  or  equal  to  one.  The  reason  is 
thought  to  be  due  to  a fluctuation  in  gas  flow  or  to  the 
electronics.  All  of  the  peak  heights  in  the  discharge 

are  corrected  for  by  K.  . 

Ar 

Ionization  Cross  Section 

From  equation  (2.16),  it  can  be  seen  that  the  ioniza- 
tion cross  section  of  each  molecular  species  in  the  source 
must  be  known.  This  total  cross  section  corresponds  to 
the  particular  energy  setting  of  the  ionizer.  It  is  desir- 
able to  have  the  ionizer  energy  2-3  eV  above  the  greatest 
ionization  energy  expected  among  the  molecular  species 
in  a mixture.  This  setting  of  the  ionizer  energy  reduces 
the  possibility  of  double  ionization  and  molecular  dissocia- 
tion occurring  in  the  ionizer.  In  this  experiment,  the 
highest  ionization  energy  is  that  of  He  (24.58  eV)  . A 
typical  setting  of  the  ionizer  energy  would  then  be  27.0  eV. 
When  the  contact  potential  is  taken  into  account,  the  meter 
reading  is  27.0  eV  plus  the  contact  potential.  The  contact 
potential  is  explained  in  the  experiment  arrangement  and 
procedute  section.  It  is  found  experimentally  that  the 
best  signal  results  are  achieved  with  an  energy  between 
65  and  70  eV. 
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The  Reduction  Process 

If  the  ionization  cross  section  for  an  additional 
species  is  known  and  the  signal  is  measured,  equation  (2.16) 


can  be  solved  for  n giving  the  total  number  density  of 

5 

species  s that  is  produced  in  the  discharge: 


_ - <PH>  o 

ns  s 

s a F 
s s 


(A.  2) 


When  the  argon  correction  factor  is  included,  the  reduced 
percentage  of  molecular  species  s becomes: 


ns  = 


(PH) sKAr 
°sPs 


(A.  3) 


Appendix  B 

Residence  Time  and  Number  Density  Calculation 


Consider  the  discharge  tube  presented  as  G in  Fig.  1. 
Assuming  Poiseuille  flow  through  this  tube,  the  radial 
velocity  profile  is  parabolic.  As  a result,  the  gas  near  the 
walls  have  a slower  velocity  and  thus  stay  in  the  discharge 
longer.  It  is  the  gas  near  the  walls  that  the  spectrometer 
samples  because  the  sampling  orifice  is  located  at  the 
wall.  In  the  following  discussion,  it  is  assumed  that 
all  the  molecules  travel  through  the  tube  at  the  same  rate. 

The  volume  of  gas  that  passes  a specified  point  in 
time  t is  given  by: 

V = vtA  (B.l) 

where  v is  the  velocity  of  the  flow,  and  A is  the  cross 
sectional  area  of  tube  (which  is  assumed  constant) . Solving 
for  the  flow  velocity  in  equation  (B.l),  the  following 
result  is  obtained: 

v - u (B-2» 

V 

— is  the  volume  flow  rate  past  a point  within  the  tube. 

This  flow  rate  is  usually  given  in  standard  cubic  centi- 
meters per  minute  (SCCM) . It  is  the  volume  occupied  by 
a given  gas  at  standard  temperature  and  pressure  (STP) . 
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Standard  temperature  and  pressure  are  defined  as  20 °C  and 


[PT 


760  mm  of  Hg  (torr)  respectively.  Assuming  an  ideal  gas, 
the  number  density  of  molecules  in  one  standard  cubic  centi- 
meter is  given  by: 

P 

_ _ o 

o kT 

o 

nQ  = 2.505  x 1019  cm-3  (B.3) 

where  P and  T are  the  standard  values  of  pressure  and 
o o 

temperature  respectively,  and  k is  Boltzmann's  constant. 

If  the  temperature  is  constant,  then  the  volume  occu- 
pied at  any  arbitrary  pressure  P is  given  by 

P 

V = VQ  (B.4) 

where  VQ  is  the  volume  occupied  at  standard  pressure.  Sub- 
stituting V from  equation  (B.4)  into  equation  (B.2),  the 
flow  velocity  becomes 


v 


P 

o 

P 


(B.  5) 


V 

where  H = — . H is  the  standard  flow  rate  given  in  SCCM. 

As  stated  in  the  text,  the  residence  time  is  defined 
as  the  length  of  time  that  a portion  of  the  gas  stays  in 
the  discharge  before  it  is  sampled  by  the  spectrometers. 
The  residence  time  is  given  by 


T - ^ (B.6) 


i 

i 

i 
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i 


where  1 is  the  length  from  the  cathode  to  the  sampling 
orifice.  Substituting  the  flow  velocity  from  equation 
(B.5)  into  equation  (B.6) , the  residence  time  becomes: 


T 


A1P 

HPo 


(B.  7) 


When  the  flow  of  gas  through  the  discharge  tube  reaches 
steady  state,  the  pressure  remains  constant  and  the  flow 
in  is  equal  to  the  flow  out.  Under  this  condition  of  steady 
state,  the  total  number  of  molecules  in  the  experimental 
region  of  interest  (i.e.,  from  the  cathode  to  the  sampling 
orifice)  is  given  by 


N = n HT  (B. 8) 

o 

where  T is  the  residence  time  of  the  gas,  H is  the  total 
flow  rate,  and  nQ  is  the  number  density  of  molecules  in 
a standard  cubic  centimeter.  If  the  flow  contains  a mixture 
of  gases,  then  H is  just  the  sum  of  the  individual  flow 
rates  of  all  the  gases  that  constitute  the  mixture: . 

i 

H = £H  (B.  9) 

s S 

Substituting  the  total  flow  rate  from  equation  (B.9)  into 
equation  (B.8),  the  total  number  in  the  experimental  region 

ts  fivan  by 
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N - n0T^T- 


(B.10) 


Using  equation  (B.8),  the  number  density  is 


n - A - n T Y 
A1  © jz  JU 

s 

n ■ no  sir  E/  Hs 


n = n P 

05i 


(B. 11) 


where  P is  the  total  pressure  of  the  gas  mixture.  The 
total  number  of  a single  component  is  given  by  nQTHg,  and 
the  number  density  by  ncTHg/Al  so  that  the  fractional  amount 
of  species  s is  given  by 


n E 


(B.12) 


Expression  (B.12)  uses  the  fact  that  the  residence  times 
of  all  the  gas  components  are  the  same. 
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Appendix  C 

Important  Reactions  and  Rate  Constants 

C02/Ar 


Reaction 

Rate  Constant  (300°K) 

Reference 

1. 

C02  + e”  CO  + 0 + 

e 

_ . . -11  3 -1 

5.4  x 10  cm  s 

* 

2. 

0 + 0 + M**  + 02  + M 

_ . , -33  6 -1 

7.2  x 10  cm  s 

6 

3. 

C2  + e"  0 + 0 + e“ 

, , -9  3 -1 

1 x 10  cm  s 

18 

C02/N2/Ar 

1. 

COz  + e”  + CO  + 0 + 

e 

, - ..-10  3 -1 

1.7  x 10  cm  s 

* 

2. 

o + o + m-»-o2  + m 

, _ . -33  6 -1 

7.2  x 10  cm  s 

6 

3. 

02  + e”  ->  O + O + e~ 

1 x lO-9  cm3  sS_1 

18 

4. 

N2  + e~  -*■  N + N + e" 

-9  3 -1 

1 x 10  cm  s 

*** 

5. 

n+n  + m-*n2+m 

„ . , -33  6 -1 

4.4  x 10  cm  s 

6 

6. 

N + 02  -*■  NO  + 0 

o c m-17  3 -1 

8.5  x 10  cm  s 

6 

7. 

N+O  + M + NO  + M 

, . -32  6 -1 

1 x 10  cm  s 

6 

8. 

N + NO  -*■  N2  +0 

, _ ..-11  3 -1 

2.7  x 10  cm  s 

6 

9. 

0 + N02  02  + NO 

. . , -12  3 -1 

9.2  x 10  cm  s 

6 

COj/He/Ar 

| 

1. 

C02  + e”  -►  CO  + 0 + 

e 

. . , -10  3 -1 

2.2  x 10  cm  s 

* 

2. 

o + o + m->-o2  + m 

7.2  x 10  33  cm**  s * 

6 

3. 

02  + e"  -*■  0 + 0 + e" 

, -9  3 -1 

1 x 10  cm  s 

18 

C02/N2/He/Ar 

1. 

C02  + e~  CO  + 0 + 

e 

„ , , -10  3 -1 

2.3  x 10  cm  s 

8 

2. 

o + o + m-*-o2+m 

_ . ..-33  6 -1 

7.2  x 10  cm  s 

6 

3. 

02  + e"  -*•  0 

+ 0 + e" 

1 X 

, -9  3 -1 

10  cm  s 

18 

4. 

N2  + e”  -*■  N 

+ N + e 

1 X 

. -9  3 -1 

10  cm  s 

*** 

5. 

N+N  + M-*- 

N2  + M 

4.4 

, -33  6 -1 

x 10  cm  s 

6 

6. 

N + 02  -*■  NO 

+ 0 

00 

• 

U! 

, -17  3 -1 

x 10  cm  s 

6 

7. 

N + O + M-*- 

NO  + M 

1 X 

. rt-32  6 -1 

10  cm  s 

6 

8. 

N + NO  ■+  N- 

+ 0 

2.7 

-11  3 -1 

x 10  cm  s 

6 

2 

♦The  average  rate  constant  calculated  by  equation  (4.8). 
♦♦All  Species. 

♦♦♦Assumed  rate  cons  taint. 


Appendix  D 

Electron  Transport  Parameters 


Drift  Velocity  ( cm/sec 


I it  I0'16  2 3 4 5 

E/N  (V  cm2) 


Figure  17.  Electron  Transport  Parameters  Used 
for  the  90%  CO2/10%  Ar  Discharge  (Ref  9) 
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Figure  19.  Electron  Transport  Parameters  Used 
for  the  10*  CO5/*0%  He/10%  Ar  Discharge 
(Ref  9) 
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